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1. Introduction 

The measurements of B c meson mass and lifetime in CDF [JTJ] and DO [§] experiments were 
the first steps in the experimental research of heavy quarkonia with open flavor. The obtained 
experimental results are in a good agreement with the theoretical predictions for the B c mass |], 

m CDF = 6.2756 ±0.0029(stat.) ±0.0025(sys.) GeV, 
m D° = 6.3000 ±0.0014(stat.) ±0.0005 (sys.) GeV, 
m£ eor = 6.25 ±0.03 GeV; 
as well as for the decay time [Q, [7|, ||, ||, [!(]] : 

t cdf = o.448lg;g||(stat.) ±0.032(sys.) ps, 
t£° = 0.475 (stat.) ±0.018(sys.) ps, 

T ±eor = ()48±005ps 

The B c meson is observed at Tevatron only in two decay modes : B c , — > J / yn and B c — > J / y + 
/I + v^j. The investigation of other decay modes is possible at at LHC, where about 10 10 events with 
B c mesons per year are expected. This huge amount of events will allow to obtain the information 
on the production cross section distributions, on the decay branching fractions, and in some cases, 
on the distributions of decay products. And absolutely new results on B c decays have already 



obtained at LHC: the B c decay to J /y + 'iil has been first observed by LHCb collaboration [11]. 

The B c systems do not have strong and electromagnetic annihilation decay modes. Due to this 
reason the exited B c systems laying below B + D threshold have the decay widths by two order of 
magnitude smaller than the widths for analogous exited states of charmonium and bottonium. All 
exited B c states after a set of radiative transitions decay into the lightest pseudoscalar state (0~). 
The lifetime of this state is comparable with the lifetimes of of B and D mesons and essentially 
differ from lifetimes of other lightest quarkonia: r\ c and rjj,. This is why B c meson provides a 
unique possibility to investigate the both strong and weak interactions. 

The main B c decay modes, such as B c — > J/y ± £v, B c — > J / y% and B c — > B, p were theo- 
retically studied in details (see, for example [|[ [12|, [T3p). For all these processes it is assumed that 
the factorization approach is valid: the decay B c — > heavy hadron + W* is followed by the decay of 
the virtual W-boson. The transition form-factors for the processes B c — > heavy hadron + W* can be 
estimated within QCD sum rules, within quark potential models, or in the framework of light-cone 
quark model. 

In our recent articles 0, [T5|] we have studied the B c decay process with several pions in the 
final state, such as B c — > J/y + nn and B c — > Bg +nn with 1 < n < 4. The factorization approach 
have been used in these calculations. The characteristics of the virtual W -boson decay have been 
adopted from the experimental data on z — y v% ± n% decays for n < 3 and from the experimental 
data on 71 mesons production in the process e + e~ — > \%. 

The modules for calculation of these decay amplitude are implemented into EvtGen program 



package [16]. Therefore the detailed simulation of these decays in the kinematical condition of 
real experiments is possible now. In this paper we describe the installation and using the developed 
modules. 
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2. Calculation technique 

To calculate the decay amplitudes 



we assume that the discussed processes 
can be represented as the decay B c — > 
heavy hadron + W* followed by the decay 
of the virtual W -boson. An additional as- 
sumption should be made that the only one of 
the constituent quarks decays weakly, mean- 
while the other quark remains the same (a 




typical diagram is presented in Fig. |l|). c 

Within this approach the amplitude of 
the process can be written in the form 

Figure 1: Typical diagram forZ? c - — > V(P) +nn decay. 

^[B c ^J/ ¥ +nn] = ^t ai (n R )J^ /v 4, 

(2.1) 

where Ew is the polarization vector of W* ; 

jifJ/v i s the B c -> J/y + W* transition vertex: 

JT M = (j/y |c 7m (1 - y 5 )b\ B c ) =%-^. (2.2) 

Vector and axial currents are equal to 

^ = (j/y\c^b\B c ) = iE» va hy Pa qpF v (q 2 ) , (2.3) 



s/n = (j/\i/\cY»Y5b\Bc) = 

e$F£ (q 2 ) +p, {EVp Bc )F* (q 2 ) +q, (e^p Bc )F A (q 2 ) , (2.4) 

where p Bc and pj/y are the momenta of B c - and //i//-mesons; 
q = pb c — Pj/y is the momentum of virtual W -boson; 

P = Pb c +Pj/ ¥ ; 

£ J Jy is the polarization vector of J/y meson; 

and F v (q 2 ), F$(q 2 ), F*(q 2 ), F*{q 2 ) and F v (q 2 ) are form-factors of B c ->■ J/y + W* decays. 

In the tree approximation the parameter ai(jUs) is equal to unity. Higher-order corrections 
lead to dependence of this factor on renormalization scale }Xr [jH|]. Numerical values for «i (jU/?) 
at different scale are calculated in [ |12| ] . For the process B c — > J/y + nn the value of jX B has been 
chosen to be equal to the mass value of the decayed ft-quark: 

a\ (mi 7 ) = 1.4 

The form-factors were calculated within different nonperturbative approaches: QCD sum rules 
jl2|], potential quark models and Light-Front quark models [JI^, |T^]. In our article we use 
exponential parametrization of these form-factors (see Tab. [[]): 

Fi{q 2 ) = F^expj-CK? 2 -^ 4 } . (2.5) 
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The decay B c — > B* + nil is described by analogy with 
the decay B c — > J/y + n7l. The same formula for the am- 
plitude is used with the other parameter set (see Tab. |l|). 

The amplitude for the decay B c — > B s + nil can be writ- 
ten in more simple form 



&/[B c ->B s + nlt] 



^ai(M*).<'4> (2-6) 



where 



Jt$> 



(B s \C (1 - y 5 ) fc| B c ) = /+ (tf 2 )^ + /_ (^r 2 )^. 



(2.7) 

The detailed information about W* decay is not needed 
to obtain the integrated branching fractions, as well as the 
branching fraction distributions on q 2 . Let us consider the 
decay B c — > J/yW* — > J /ynft as an example: 



B c — >J/y + nn 




W) 


c\ 




Ao 


5.9 


0.049 


0.0015 


A + 


-0.074 


0.049 


0.0015 


V 


0.11 


0.049 


0.0015 




B c ^ 


B* + nll 




W) 


c\ 


C2 


Ao 


8.1 


0.30 


0.069 


A+ 


0.15 


0.30 


0.069 


V 


1.08 


0.30 


0.069 




B c ^ 


B s + nK 






Mo) 


c\ 


C2 


U 


1.3 


0.30 


0.069 



Table 1: Form-factor parameters for 
different SR form-factor sets. 



dT(B c ^J/y@) 



1 



2M 2 

where Lorentz-invariant phase space is defined according to 

rfO (G -> pi ...p n ) = (2x) 4 8 4 (Q - £ Pi ) n 
Using the following recurrent expression for the phase space 



e v v d^(B c —?-J/ynx) . 



d 3 pi 



2Ei{2%) 



da 2 

d&(B c -tJ/vtW* -^J/ynn) = -^-d<S>(B c -> J/yW*)d<Z>(W* 

271 



nil) 



one can perform the integration over phase space of the final state n%: 

J^Jd® (W* -> Ml) eje™* = (q^qy - q 2 g^ v ) p T {q 2 ) + q^vpL (<? 2 ) , 



(2.8) 



(2.9) 



(2.10) 



(2.11) 



where spectral functions pji (q 2 ) are universal and can be determined from theoretical and ex- 
perimental analysis of some other processes, for example T — > v t n% decay or electron-positron 
annihilation e + e~ — > Ml. (See Jl4], p3j, E0] for details.) It is worth to mention, that due to the 



vector current conservation and the partial axial current conservation spectral function pi for n > 2 
is negligible in almost whole kinematical region, so it can be neglected in the estimations for n>2. 
For the purposes of our article, however, a more detailed description is required for the multipion 
final state. 

In the framework of resonance model the decays of W* -boson can be described in terms of 
virtual p- and ai-mesons exchange (see typical diagrams presented in Fig. |2|). 

If there is one 7i-meson in the final state (see Fig. ^|a), the vertex of W* — > 71 transition can be 
written in the form 

(7i + \(d~u) v ^ A \W)=f K k^ (2.12) 
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Figure 2: Typical diagrams for W* — > nn decay in resonance approximation. 



where k is the momentum of 7T-meson and f % rj 140 MeV is its coupling constant. In accordance 
with this interaction vertex the effective polarization vector ef[ in this case has the form 

= fAn/ml (2.13) 

Note, that this vertex violates the axial current. 

The decays of virtual W -boson into multipion final states are described within resonance model 
in terms of virtual p- and ai-mesons exchange (see typical diagrams presented in Fig. $p,c). 

The W* — > 7T + 7r° decay is saturated mainly by contributions of virtual p- and p'-mesons (see 
Fig. |||b). The corresponding effective polarization vector can be written as 

ef =F p (q 2 )(k 1 -k 2 )^ (2.14) 



where are 7C-mesons momenta and F p {q 2 ) is the p-meson form-factor (see [21 ]). The difference 
in 7T°- and 7T + -meson masses is neglected, thus the virtual W boson in this decay has a transverse 
polarization. It should be noted that the width of the p meson must be taken into account. 

The W* — > 37T-transition is saturated mainly by W* — > a\ — > p% — > 3n decay chain. Following 



[ ]21[ ] one can write the effective polarization vertex in this case as 



e f = -*Tr F «(<? 2 ) { b p( s 2)Vi^ +B p ( Sl )V 2 ^}, (2.15) 

where 

Vi2n=k l ^-k 3 ^-q^ l (2.16) 

and 

sia = (l-h,2) 2 - (2.17) 



Parametrization of B p (s) function is presented in [plQ. It can be clearly seen, that if one neglects the 
difference between charged and neutral 7C-meson masses the above expression in transverse and the 
axial current is conserved. In EvtGen package this transition was realized already in TAUHADNU 
model. 



5 



B c -> J/\j/(B s ,B*) + nit Decays 



Alexander Berezhnoy 



3. Software structure 



EvtDecayBase 



EvtDecayAmp 



EvtBcNPi 

EvlBcToNPi.cpp 
EvtBcToNPi.hh 

public: 

string getNameO; //{return "BC_NPI";} 
void init(); 
void initProbMaxf); 
Void decay(EvtParticle *); 
private: 

EvtComplex Fpi(ql,q2); 
double pi3G(double Q2); 



For the Monte-Carlo simula- 
tion of the discussed decays the 
generator EvtGen is used, which 
is the part of LHCb software envi- 
ronment GAUSS. EvtGen is writ- 
ten completely on C++. In order to 
add a new decay to it one should 
just create a new class, that de- 
scribes this decay. The base class 
is EvtDecayAmp, where proto- 
types of all necessary for the model 
functions are given. The most 
significant of these functions are 
init(), initProbMax ( ) and 
decay () (see Fig. ||]). 

The method init() per- 
forms the initialization of the de- 
cay model and reads its parameters. The necessary parameters are stored in the so-called .dec-file 
(the exact name of this file can be determined by the user). The detailed description of the .dec-file 
format can be found in the EvtGen documentation [16]. 

As it is seen from Fig. |] and || the results of Monte-Carlo simulation for the decays B c — > 
J/\jf + 2n and B c — > J/y + 3n are in good agreement with the theoretical predictions [|l4]]. 



EvtBcPsiNpi 



EvtBcBsStarNpi 



EvtBcBsNpi 



Figure 3: The program structure. 



d Br 1<T 2 
d i/ 1 GeV 2 




0.5 1.0 1.5 



Figure 4: The distribution over the squared in- 
variant mass q 2 = m^ n in B c — ¥ J /y + 2n decay 
generated within BC_NPI-model (histogram) in 
comparison with the predictions Jl4|]. 



f/Bi 1 




Figure 5: The distribution over the squared in- 
variant mass q 2 — in B t — s- J /y + 37T decay 
generated within BC_NPI-model (histogram) in 
comparison with the predictions [jjjj]. 
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4. Results 

• The branching ratios for the decays B c ^-J/ y{B s ,B*) + nn have been estimated in the frame- 
work of the factorization approach. 

• The width of intermediate mesons must be taken into account. 

• The package for simulation of these decays within the standard LHCb software GAUSS has 
been developed. 

• The experimental LHCb data on B c — > J /y + 2>n decay JTT| ] can be described satisfactorily 
within the model. 
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